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Recently it was shown that quantum corrections to the Newton potential 
can explain the rotation curves in spiral galaxies without introducing the 



O ! Dark Matter halo. The unique phenomenological parameter au of the the- 

^ , ory grows with the mass of the galaxy. In order to better investigate the 

mass- dependence of au one needs to check the upper bound for au at a smaller 
^ scale. Here we perform the corresponding calculation by analyzing the dynam- 

>■ ■ ics of the Laplace-Runge-Lenz vector. The resulting limitation on quantum 



corrections is quite severe, suggesting a strong mass-dependence of au. 



Introduction 



It is a common belief nowadays that the General Relativity (GR) is not the ultimate 
^ ! theory of gravity. One reason for this is that the relevant solutions of GR, such as the 
spherically symmetric solution and the homogeneous and isotropic one, both manifest 
singular behavior in their extremes. In the first case the space-time singularity is in the 
center of the black hole and in the second case it is in the initial instant of the universe "Big 
Bang" . In both cases the singularity is surrounded by a very small space-time region with 
very high magnitudes of curvature tensor components. This makes perfectly possible that 
the higher derivative terms in the gravitational action may change the geometry in such 
a way that the singularities should disappear. The importance of higher derivative terms 
in the gravitational action is due to the fact that they are requested for constructing 
a renormalizable theory of matter (including Standard Model) on curved background 
(see, e.g., [1] and further references therein). The effects of higher derivative terms on 
singularities were discussed in the cosmological (see, e.g., |2l[3l|4]) and black hole (see, e.g.. 
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[5l|6]) settings and there are serious reasons to consider the possibihty of erased singularity 
due to the higher derivative terms in the classical action and quantum corrections |7]. 

Despite the discussion of higher derivative terms and their possible effect on singu- 
larities is interesting, we are much more curious about possible modifications of gravity 
and, especially, about possible quantum effects at low energy scale, where the observa- 
tions are much more real. At low energies the effect of higher derivative terms is usually 
assumed to be Planck-suppressed and one has to deal with, e.g., quantum corrections 
to the Hilbert-Einstein action. For the higher derivative section of the theory, direct 
calculations of the low-energy quantum contributions were recently performed in [SI [9] 
using the approximation of linearized metric on flat background. At the same time it 
was demonstrated in [8] that this kind of calculational technique is useless for deriving 
quantum corrections for the cosmological and Hilbert-Einstein terms. Therefore, despite 
this part is much more interesting from the physical point of view, here we have much less 
achievements. As far as the subject looks quite relevant for applications, it is worthwhile 
to try some phenomenological approach in this case. 

In what follows we will consider the low-energy quantum corrections to the Newton 
potential using the renormalization group technique and the identification of renormal- 
ization scale which was recently proposed in [10] on the phenomenological basis and then 
justified theoretically in [TT]. This identification of scale includes some uncertainty which 
is measured by a dimensionless phenomenological parameter au. It was shown in p2] 
that values of au of the order 10~^ can provide the detailed and precise explanation for 
the rotation curves for a small but quite representative sample of spiral galaxies. It is 
remarkable that the mentioned parameter is steadily growing with the increase of the 
mass of the galaxy. Therefore, it is natural to expect a much smaller value for au for 
much smaller astrophysical systems, such as starts. Our purpose is to set an upper bound 
on au for the Solar system, by using a very efficient approach based on the dynamics 
of the Laplace-Runge-Lentz vector (see, for instance, [12] and references therein for the 
introduction, for an interesting historical review, see [T3]). 

The paper is organized as follows. In the next section we present very general argu- 
ments about unique possible form of renormalization group running of Newton constant 
G in the IR region. In Sect. 3 a brief introduction and necessary information about the 
Laplace-Runge-Lentz vector in the almost Keplerian problem are presented. The numer- 
ical estimates and the upper bound for au from the Mercury precession are derived in 
Sect. 4. In Sect. 5 we draw our conclusions. 
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2 Quantum effects and Newton potential 



One of the most powerful techniques for evaluating quantum corrections is renormalization 
group. So, let us check out what is the renormalization group equation in the low-energy 
gravitational sector. The unique relevant parameter at the astrophysical scale is the 
Newton constant G, and at the quantum level it becomes a running parameter G(/i), 
where /i defines a scale. The problem of identifying fi with some physical quantity will 
be discussed later on, and now we concentrate on the dependence which is always 

governed by the corresponding renormalization group equation ^{dG/dfi) = Pg- 

Consider an arbitrary quantum theory with gravity. It can be, for instance, some 
quantum theory of the gravitational field or quantum theory of matter fields. Every kind 
of quantum theory can be characterized by the massive parameters m which define scale. 
There may be, of course, more than one such parameter, so let us consider, for the sake 
of generality, the whole set {rrti}. For example, the elements of the set {mj} can be the 
masses of all particles or fields which are present in the given quantum theory. 

Let us present general arguments about the possible form of the running Newton 
constant, G{fi). Using dimensional arguments we can establish the unique possible form 
of the renormalization group equation, 

/x^— = J2 AjmiTRj = 2pMI, G-\fio) = Go' = Mj, . (1) 

particles 

In particular, in the SM-like or GUT-like theory, at the one-loop level, one has 

particles fermions scalars 

where the fermion masses were denoted by m/ and C,s is the nonminimal parameter for 
the scalar with the mass m^. 

At the same time, it is important to stress that Eq. ([1]) is not just a one-loop equation, 
but it is valid at any loop order and in any theory which is capable to produce the renor- 
malization group equation for G. In general, beyond the one-loop level, the coefficients 
Aij in Eq. ([1]) depend on coupling constants which are present in the theory. 

One can rewrite Eq. ([1]) as 

/.^^^ = -2.(G/Go)^ (3) 

Solving this equation we obtain the universal form of the scale dependence for the Newton 
constant. 
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Here the word universal means that either G is not running (and this means there 
are no quantum effects in the low-energy gravity sector) or such running is given by Eq. 
(jlj). From this perspective it is not a surprise that Eqs. (jS]) and (jlj) can be met in such 
quahtatively different approaches as higher derivative quantum gravity [HI [151 [Ml [IZ] , 
quantum theory of matter fields on curved background [HI H] and quantum theory of 
conformal factor [19]. The same equation shows up also in the phenomenological 
approach based on the hypothesis of the Appelquist and Carazzone - like decoupling for 
the cosmological constant and conservation law for the quantum corrected gravitational 
action [20]. The reason behind these occurrences is that any other form of G{^) would 
be in conflict with very simple (and hence very safe) dimensional considerations and also 
with the covariance arguments which play a very significant role here [7]. 

The next problem is how to identify fi. As usual, this identification depends on the 
physical problem under discussion and there is no universal solution. In the case of the 
gravitational field of a point-like mass, the most natural choice is fi ~ 1/r, where r 
is the distance from the mass position. This choice of the scale identification has been 
used in various publications [211 ESI ESI EO]- In particular, this identification enables one 
to roughly explain the fiat rotation curves for the point-like model of the spiral galaxies 
[211 l22l [20] by directly using Eq. (jl]). Moreover, one can achieve very good and detailed 
description for the good sample of rotation curves by using Eq. (jl]) together with more 
sophisticated identification of the scale. 



where the value of $0 is irrelevant and $Newt is the Newtonian potential computed with 
the boundary condition of it being zero at infinity. In Eq. a is a phenomenological 
parameter which should be defined from fitting to the observational data. Let us note that 
the same dependence ([5]) can be also obtained from the regular scale-setting procedure 
[24], which was applied to the present case in [TT] . 

It is easy to see from Eqs. (jl]) and ([5]) that a always shows up as a factor in the 
product au. It turns out that the fit with the observational data is perfect (for a sample 
of nine galaxies) if we assume that the product au is about 10"'' and, moreover, this 
product grows up with the increasing of the mass of the galaxy [10]. Indeed, this is a 
very nice feature, because then one may hope that the effect of "corrected" Newton law 
would be very weak at the scale of the Solar system, which has a mass of many orders of 
magnitude smaller than the one of a galaxy. The purpose of the present paper is to make 
the last statement quantitative, that is to set an upper bound on the value of au inside 
the Solar system. The best available data here are about the precession of the perihelion 
of Mercury, so it is sufficient to deal with these data only. 




(5) 
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3 Laplace-Runge-Lenz vector 



The method that we will use to calculate the precession in the orbit of the Mercury, due 
to quantum effect in the Newtonian potential, is based on the known Laplace-Runge-Lenz 
vector. Therefore we shall do a discuss briefly about it and we will see as we can use it 
to calculate the precession velocity. 

Consider first the motion of a particle of mass m for the non-perturbed Kepler's 
problem, when there is a single Newtonian force acting on the particle, 

k 

FjVe«;t = -^f, (6) 

where k — GoMm. The Laplace-Runge-Lenz vector (LRL) is defined as 

A ^ p X £ - mkr , (7) 

where £ = r x p is the angular momentum vector. One can show that in the case of the 
non-perturbed Kepler's problem LRL is a constant of motion, namely, 

dA dp „ d£ , dr ^ 

— ^ ^ X £ + px — -mk— ^ 0. (8) 
dt dt dt dt ^ ^ 

Furthermore, LRL has some important relationships with other constants of motion. For 
example, £-A — 0, hence LRL always remains in the plane of the orbit. Another important 
relation concerns the total energy of the particle E and to its angular momentum 

(2E£^\ 
1 + . (9) 

remembering that for the Kepler problem the eccentricity of the orbit is related to the 
energy and angular momentum as 



so that we can write the modulo of the LRL vector as 

|A| = mke . (11) 

One can see that the magnitude of the LRL vector measures the eccentricity of the orbit. 
Moreover, one can arrive at the equation of the orbit by taking scalar product of the LRL 
with the position vector, namely 

r| A| cos((^ — (^o) — {px £ — mkr^ — ^ — mkr , (12) 
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where (po is the angle between A and polar axis. After some simple manipulations we get 

£'^/mk 



(13) 



One can see that A is pointing to the direction of symmetry of the orbit, which can be 
limited or unlimited. Then it is convenient to choose the polar axis in the direction of the 
LRL vector A, e.g., by choosing ipo = 0. Then A = mks'k and 



e/mk a(l 



e 



2^ 



e COSV9 L + e cosip 



(14) 



where a is the major semi-axis of the ellipse. 

Once we know the properties and the interpretation of the LRL vector, we are able 
to deal with the Kepler's problem when a small perturbation is introduced. In this case 
the new orbit will be very similar to the old one, however there will be a precession. In 
other words, the particle has an approximately elliptical orbit, but in such a way that the 
major semi-axis slowly rotates. The velocity of this rotation is called precession velocity. 
A comparison between theoretical predictions and available experimental (observational) 
data may give valuable information such as, for instance, upper bounds on relevant pa- 
rameters which are present in the perturbing force. In our case, we shall follow this line 
by using the well-known data for the Mercury precession. 

Consider a particle of mass m which moves under the action of a total force 

F = -Af + f. (15) 

Here f is a small perturbation force (|f | ^ kjr^\ which can be non-central, in principle. 
Some formulas below are general, but in fact we are mainly interested in the central 
perturbed force f = — gradM(r), where the perturbing term for the potential energy can 
be derived from the Eqs. ^ and (151) to be pUf [TU] 

$ = $Newt + u{r) , u{r) = — - — , (16) 

where SG = G{fi) — Go- 
Let us present, for the sake of completeness, the short review of the derivation of Eq. 
( fT6l) . One can start from the expression for the action with variable G, 



— / a xJ—q ——^ 



where /i is supposed to depend on some energy features of the gravitational field. We 
assume that G has a very weak deviation from the constant value Go = ^/Mpi, namely 
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G = Go{l + k). In what follows we will consider k to be a very small quantity and hence 
keep only first order terms in this parameter. Our purpose is to link the action ( JT71) with 
the usual one of GR, with Go instead of G{fi). For this end we perform the conformal 
transformation [20j according to 

9^,u = ^ g^^u = (l- K)g^^ . (18) 

The derivatives of k emerge in the transformed action, but only in second power, hence 
they may be neglected. Then, in the linear order in k, the metric cjf^i, satisfies Einstein 
equations with constant Go, and the nonrelativistic limit of the two metrics gfj_u and g^i, 
is related as 

gQQ = -1- — , hence goo = -1 — . (19) 

Here $Newt is the usual Newton potential and $ is an apparent potential corresponding 
to the nonrelativistic solution of the modifies gravitational theory (fT7|l . Due to the (flQll . 
we arrive at [201 [IH] 

= $Newt + Trl^ = "^Newt + TTT^ > (20) 
Z / (jro 

which is nothing else but f|T6|) . Now we can come back to the analysis of this formula. For 
the Solar system case we can use the relation $Newt = —k/r. Hence, the identification 
dSJ effectively coincides with the ~ 1/r one of the Refs. [HI ESI ESI ED] ■ 

One can easily prove the following statement: For the case of a central perturbation 
force f = —grad 'u(r), the magnitude of the LRL vector varies according to the relation 

dA^ = -2m f du{r) = -2m u (r) dr . (21) 

The last relation shows that the magnitude of the vector A varies as the distance r 
varies, but in such a way that it takes equal values for equal values of r. If we restrict 
the discussion by quase-elliptic orbits, which are restricted in space, we have ri < r < r2, 
where ri and r2 are the corresponding turning points, then the modulus of the Laplace- 
Runge-Lenz vector A will assume the same value whenever the distance r is the same. 
Hence, the dynamics of A is perfectly useful for evaluating the precession of the orbit 
in the quasi- Newtonian case. The variation of the magnitude |A| is of the first order 
of magnitude in the small perturbation u{r). Therefore, we can completely neglect this 
variation when evaluating the precession of the LRL vector, because this precession is 
also of the first order in u{r). In what follows we will assume, for the sake of simplicity, 
that |A| is constant even when the small perturbations are present. 
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The time variation rate of the LRL vector is given by 



dA 



i X i + px(rxf). 



(22) 



As we have seen above, the LRL always points towards the symmetry axis of the orbit. 
In the Kepler problem, for elliptical case, this symmetry axis is the major semi-axis of 
the orbit. Therefore, we can calculate the velocity of precession of the orbit by simply 
computing the velocity of precession of the LRL vector. In the present work we are 
concerned with the quantum correction to the Newton gravitational potential. Such 
correction is expected to be a very small quantity, hence it is completely fair to employ 
a perturbative approach. With this consideration in mind, we compute the time average 
of the velocity of precession of the LRL vector for one period of the unperturbed orbit. 
One can recall that both £ and A are constants of motion for the unperturbed orbit. 
Therefore, for the sake of our calculation we can simply use £ and A, instead of {£) and 
(A). Taking, then, the time average of the time derivative of the LRL vector over the 
period of the unperturbed orbit, we obtain 



where A is the LRL vector in the unperturbed case, and we denoted by Q the time average 
value of the velocity of precession of the LRL vector. Let us note that, in the l.h.s. of the 
Eq. fl23l) . there is a quantity which is of the first order in the perturbing force. That is 
why the sign of averaging can not be omitted here. 

In Eq. ( l23l) the time-averaging of a function F means 



where r is the period of the non-perturbed motion (for the closed orbit case. However, 
since we are mainly interested in the trajectory, r{ip), it is convenient to perform a change 
of variables and trade the time integration for the angular integration, namely. 



We note that for a central perturbative force f , the angular momentum £ is constant 
and the second term in the r.h.s. of Eq. fl22|) vanishes. Then one can write, in the 
Cartesian basis. 




(23) 




(24) 




(25) 




= (f X €) = (/(r) cosy?) X X + {f{r) simp) 



yx£. 



(26) 
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Since / depends only on the distance r and the non-perturbed orbit is symmetric with 
respect to the major semi-axis, the second term in the r.h.s. of ( l26l) vanishes and we 
obtain 

dA\ ^ . , ^ ( f(r) cos cp) , , , 

— )=^^xA, where n = -^l^^ — 27 
at I mke 

describing the average precession of the orbit. 

4 Solar system tests for logarithmic term 

Let us now apply the method developed in the previous section to the case where the 
disturbing force is due to quantum effects. As we have already seen in Sect. 2, the 
corresponding additional term for the gravitational potential is given by 

, , m(? 5G , , 

«W = ^^, (28) 

where, according to Eqs. (j4]) and ([5]), we meet logarithmic dependence 

G($Newt) = — • (29) 

1 + z/a In ($Ncwt/$o) 
From the previous equation, the perturbing force acting on the particle is given by 

f = - ^ gradG($Newt) = - -— • (30) 

From Eqs. ( 127|) and ( 130|) . we obtain the following average precession rate. 



vac t I cos if 
ke 



l + ^«log (^) 



-2\ 



rke 1 Jo [l + 2ua ln(ro/r)]' ' 

One can remember that already at the typical galaxy scale ua oc 10"'' and that the 
expected bound for the Solar system should be essentially smaller. Therefore it is justified 
to use an approximation z/a ^ 1 and keep only first order terms in this parameter. Then, 
after some simple calculations we arrive at the following expression for the absolute value 
of the precession velocity: 



iZ ~ ; / r(09)COS09 

rke In 



1 - 2ua In 



dip 
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ua mc ' 







r{ip) cos {pd(p 

"iR in 

dip . (32) 



rke 

vamc^a{l — e^) f'^'^ cos 99 



rke Jq 1 + £003^9 

The integral in the above expression can be easily solved and as a result we find 



_ 27rac^ au[Vl^-{l-e^)] 

~ • ^^''^ 

One can use Eq. ( l32l ) for deriving an upper bound for the parameter av in the Solar 
system. For the case of the precession of Mercury, the uncertainty in the measurement of 
the velocity of precession is 0.45" per century. 

After some simple calculations we arrive at the upper bound, 

z/a<10~^^ (34) 

where we have used the following values: 

G = 6.67 ■ 10'^^ Nm^kg-^ , M = 1.98 ■ 10^° kg , c = 3 ■ 10^ ms'^ 
a = 6.97 ■ 10^°m, r = 0.241 years and e = 0.2056 . (35) 

Here r, a, e are the period of rotation of Mercury, the major semi-axis and eccentricity of 
its orbit. 



5 Conclusions. 

We have considered the running of the Newton constant G in the framework of a general 
Renormalization Group approach. The covariance and dimensional arguments lead to a 
unique possible form of such running, which can take place in all loop orders. The beta- 
function for G has one arbitrary parameter u which depends on the details of the given 
theory. Vanishing u means there is no running at all, this means there are no relevant 
quantum corrections into the low-energy sector (Hilbert-Einstein) of the gravitational 
action. 

Assuming that u is non-zero, one can try to derive upper bound for u from different 
gravitational observations. In the recent paper [10] it was shown that the identification 
of the renormalization group scale fi ~ ($Newt)" provides an excellent fit for the rotation 
curves of the galaxies, with the product au being about 10~^ and moreover steadily 
growing with the increase of the mass of the galaxy under consideration. No Dark Matter 
is requested for this fit. In the present work we have derived an upper bound for au in the 
Solar system using the method based on the dynamics of the LRL vector. The maximal 
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possible value we have obtained for ai/ is 10~^ , which implies qualitatively the same form 
of a running that was predicted in pTj on theoretical basis. 

Indeed, both values au oc 10"^ and 10"^'' have to be seen as maximal ones, representing 
upper bounds. In case of rotation curves one can admit certain amount of a Dark Matter 
(with fi^)^/ smaller than usual), which should be helpful to explain other observations 
(LSS, CMB, BAO etc). Let us note that the first paper exploring the possibility of an 
alternative concordance model with taking into account quantum corrections is under 
preparation [25]. In the case of a quasi- Newtonian potential in the Solar system we 
can also see the value 10~^^ as an upper bound, such that the real value of au can 
be much smaller than that. However, it is definitely remarkable that the two different 
observations produced results which are consistent with each other and also with the 
theoretical prediction of [llj. 
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